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Abstrac  t 


Two-dimensional  inviscid,  unsteady  transonic  flows  with  shock  waves 
in  asymmetric  channels  with  arbitrary  wall  shapes  are  investigated,  using 
the  method  of  matched  asymptotic  expansions.  With  the  exception  of  thin 
regions  in  the  neighborhood  of  the  channel  throat  and  the  neighborhood  of 
the  shock  wave,  solutions  found  using  linearized  governing  equations  are 
valid.  In  the  regions  enclosing  the  shock  wave,  an  inner  solution  which 
satisfies  the  shock  jump  conditions  and  matches  with  the  outer  solutions, 
is  presented.  The  shock  shape  is  found  as  part  of  the  solution,  which  is 
obtained  numerically  using  the  method  of  integral  relations.  A composite 
solution,  uniformly  valid  throughout  the  channel,  and  the  relation  between 
the  instantaneous  shock  wave  position  and  back  pressure  far  downstream 
are  presented. 


INTRODUCTION 


The  analysis  of  two-dimensional  transonic  channel  flows  w'ith  shock 

waves  has  recently  been  the  subject  of  several  papers,  ^ ^ culminatine  in  a 

description  of  solutions  found  for  unsteady  flows  with  arbitrary  initial  con- 

f 

ditions  and  arbitrary  wall  shapes.  However,  the  solutions,  no  matter 
whether  the  channels  are  symmetric  or  asymmetric,  are  valid  for  stream- 
line curvatures  small  enough  that,  to  the  scale  of  the  channel  width,  the 
shock  waves  are  planar.  In  this  paper,  on  the  ■ ther  hand,  a lower  order 
(larger)  curvature  is  considered  in  an  a«vonmetrii-  channel,  so  that  strong 
cross  channel  gradients  in  velocity  and  p-essure  are  found  in  the  lowest 
order  solutions.  As  a result,  the  shock  wave  i.s  not  planar;  its  shape  is 
unknown  and  must  be  found  as  part  of  the  solution. 

Two  approaches  have  generally  been  used  in  constructing  analytical 
solutions  for  transonic  channel  flows  with  shock  waves.  In  one,  similarity 

solutions  of  the  transonic  sn'.all-disturbance  equation  are  sought  through 

1 ^ Z 

various  transformations,  resulting  in  solutions  for  steady  or  unsteady 
flow's.  The  solutions  thus  obtained,  however,  satisfy  only  very  special 
boundary  conditions  which  may  or  may  not  correspond  to  the  physical  prob- 
lem under  consideration,  and  the  solutions  are  valid  in  a thin  neighborhood 
of  the  sonic  line.  In  the  other  approach,  the  method  of  matched  asymptotic 
expansions  is  used  to  derive  governing  equations  for  each  order  of  approxi- 
mation from  the  general  inviscid  flow  equations,  for  various  regions  in  the 
channel.  Adamson,  Messiter  and  Richey^  showed  that  in  the  main  part  of 
the  channel,  this  method  gives  a systematic  derivation  of  the  power  series 
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expansion  solution  postulated  by  Szaniawski  and  used,  for  example,  in  th^^ 

3 

paper  by  Kopystynski  and  Szaniawski.  In  this  method,  arbitrary  wall 
j shapes  may  be  specified,  and  in  the  extension  to  unsteady  flows,  ^ arbitrary 

r 

1 

initial  conditions  may  be  considered.  It  is  this  approach  which  is  used  in 
the  present  paper. 

In  the  followimj,  the  flow  is  taken  to  be  two-dimensional,  compres- 
sible, and  transonic.  The  gas  is  assumed  to  follow  »he  perfect  gas  law  and 
to  have  constant  specific  heats.  Upstream  of  the  shock  wave  the  flow  is 
considered  tv  !)e  ir rotational,  and  the  Reynolds  number  is  taken  to  be  large 
enough  that  viscous  effects  are  negligible.  The  unsteadiness  in  the  flow  is 
assumed  to  arise  as  a result  of  disturbances  impressed  upon  the  flow  down- 
stream of  the  region  under  consideration.  A more  detailed  descriptioi  of 
the  analysis  which  follows,  may  be  found  in  Reference  7. 

1 

CHANNEL  FLOW  SOLUTIONS 

The  solution  is  written  in  terms  of  time  dependent  perturbations 

from  a steady,  uniform,  ir  rotational , sonic  flow.  The  coc<rdinate  system 

and  notation  used  are  illustrated  in  figure  1.  Velocity  components  u and  v 

and  the  sound  speed  a are  made  dimensionaless  with  respect  to  a“',  the 

sound  speed  in  the  undisturbed  flow  (overbars  indicate  dimensional 

quantities),  space  variables  x and  y are  referred  to  L,  the  channel  half 

width  at  the  minimum  area,  and  the  time  T is  made  dimensionless  with 

respect  to  L/a"'.  The  pressure,  density,  and  temperature,  denoted  by 
A 

P,  p,  and  T respectively,  are  made  dimensionless  with  respect  to  their 
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values  in  the  undisturbed  flow.  The  regicm  being  considered  is  that  for 
which  X = 0(1)  and  y = 0(1);  that  is,  it  extends  downstream  of  the  throat  for 
distances  of  the  order  of  the  channel  half  width.  We  introduce  a small 
parameter,  e <<  I , which  is  of  the  order  of  the  typical  percentage  difference 
between  the  flow  velocity  and  its  critical  sonic  value.  Thus  u-  1 -0(e)  for 
X = 0(1).  The  flow  entering  the  region  under  consideration  is  irrotational 
and  boundary  layer  effects  are  negligible.  Moreover,  because  the  flow  is 
transonic,  shock  waves  are  weak  enough  that  it  can  be  shown  that  one  can 
define  a velocity  potential  to  the  order  considered  here. 

Unsteady  flows  may  be  characterized  hy  prescribing  the  relative 
orders  of  the  characteristic  time  associated  with  the  imposed  flov.  distur- 
bance, T , and  the  characteristic  flow  time,  L/a*.  In  ttiis  present  work, 
c h 

the  so-called  slowly  varying  time  regime  is  coirsidered,  wliere  L/a  . 

Hence,  a parameter  t is  introduced  as  follows, 


T = T t 


T T^^/(  L/a''» 


(la,b) 


where  x » I,  and  t 0(1).  The  slowly  varying  tinie  regime  covers  a range 
of  values  for  T^j^  which  is  found  in  many  physical  problems.^  The  relation- 
ship between  x and  e depends  upon  the  relative  orders  of  T^  j^  and  the  time 
it  takes  a pulse  to  move  upstream  from  the  point  at  which  a disturbance 
originates,  to  the  location  under  consideration,  c.g.  , the  point  at  which  a 
shock  wave  is  located.  Thu.s,  a disturbance  pulse  travels  upstream  at 
sonic  velocity  relative  to  the  flow  and,  therefore,  at  an  alisolute  velocity 
Up  - 1 - u = 0(e).  Hence,  the  time  required  for  the  disturbance  to  travel 
a distance  L is  0(  L/e  iT^)  and  the  desired  time  ratio  is  ( H/e  a ' ) = x /e  . 
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It  is  seen,  therefore,  that  if  the  two  times  are  of  the  same  order,  t = 0(el. 
This  is  the  case  considered  in  reference  6.  Here  we  consider  the  case  where 
» L/  (f  a *) , and  in  fact  define  t as , 

T = (k  e ^ ' ( z 

where  k is  a constant  of  order  unity.  Physically,  this  means  that  flow 
disturbances  are  communicated  "instantaneously”  to  the  shock  wave. 


The  governing  equations  are  the  gas  dynamic  equation  and  the 
Bernoulli  equation: 


(a^ t (a^  -^TT  ' ^^x^xT  ' ^'^y^yT  = « ^^a) 


1 


<e>t  + ^y)  + - i)  = (Y  f i)/z{y  - 1) 


(3b) 


v/here  subscripts  x,  y and  T indicate  partial  differentiation;  it  can  be  shown' 
that  the  constant  on  the  right  hand  side  of  equation  (3b)  is  valid  to  and 
including  terms  of  order  , the  highest  order  terms  needed  here.  Since 
the  solutions  are  written  in  terms  of  perturbations  from  a uniform  flow,  it 
is  convenient  to  write  $ as, 

^(x,y,T)=x-t4)(x,y,T)  (4) 

where  d)  is  the  perturbation  potential.  Then,  using  equation  (3b)  for  a^  and 
equation  (4)  for  equation  (3a)  may  be  written  as  follows: 

(j,  = (Y+  l)[((j>  + (j)V2)4)  + {4>V2)4>  ] 4 2 (4>  + 4.  4>  )4> 

yy  X X XX  y yy  X X y xy 


4 (V-  1 )[(())  4 4,V2)c()  + (4.V2)4>  ] + 2 (1  +4)  )4>  „ 

X X yy  y XX  x xT 

4 (Y-  l)4>^(4>^^  + 4>yy)  + 2 4>y4>^^  4 


(^) 
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The  boundary  conditions  at  the  walls  of  the  channel  are,  for  stationary 

walls  , 

(x , y , T)  = y'  [ 1 + 4>  (x , y , T )]  (6 ) 

y w w X w 

where  y'  = dy  /dx.  For  the  asymmetric  channels  under  consideration, 
w w 

and  with  the  coordinate  system  as  showTi  in  figure  1,  the  wall  shape,  y^  is 
written  as 

y = ± 1 + € f,(x)  ±e'f2(x)  (7) 

w 

where  fj(0)  = f2(0)  = f]'(0)  = f2'(0)  = 0 and  where  fi"(x)  and  faU^)  are  continuous 
and  nonzero  at  x = 0.  The  upper  and  lower  signs  refer  to  the  upper  and  lo^er 
walls,  respectively.  It  is  seen  that  y = e fi(x)  is  the  equation  for  the  channel 
centerline  and  ± (1  +e^f2(x))  are  terms  which  describe  walls  symmetric  with 
respect  to  the  centerline,  so  that  e is  a measure  of  the  radius  of  curvature 
of  the  channel  at  x = 0,  the  channel  throat.  Thus,  the  channel  curvature  is 
an  order  of  magnitude  larger  than  that  considered  previously,  and  as  will 
be  seen,  this  leads  to  completely  different  flow  structures  and  shock  shapes. 
An  asymptotic  expansion  of  <))  is  written  as  follows: 

4>(x,  y,T)  = e4>i(x,y,t)4e^(|)2(x,y,t)-fe*4)j(x,y,t)  + (H) 

Substitution  of  equations  (1),  (2)  and  (8)  into  equation  (5)  yields  the  following 
equations  for  4>i  i 4>z>  ^.nd  <|)3  : 


syy 


yy 


4>  = [(V  + 1)(|)^  /2  + 4>M 

zyy  ix'  ^ly  X 

[(V+l)<j)  4>  + 2 (())  4>,  + k<t)  ) + ^ (V -7)  (\+ 1)4)3 

lx  2x  ly  2y  It  3 2 lx 

+ V4>  4>M  -(v  + 1)4>34>  /2 

lx  ly  X ly  ixx 


(9a) 

(9b) 


(9c) 
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The  boundary  conditions  for  equations  (9)  are  found  by  substituting  equations 
(U»  (2),  (7)  and  (8)  into  equation  (6).  Thus,  to  first  order,  it  is  found  that 

± 1 , t)  = fi'(x)  (10) 

where  the  prime  indicates  differentiation  with  respect  to  x.  The  solution 
to  equation  (9a)  subject  to  the  boundary  conditions  in  equation  (10)  is 

4>i  = fi'(x)  y + hi(x,  t)  (11) 

where  hi(x,t)  is  a function  of  integration.  If  the  solution  for  6i  , equation  (11), 
is  substituted  into  equation  (9b),  the  resulting  equation  may  be  integrated  to 
give  a solution  for  d)2.  Thus, 

4.z(x,y.t)  = fi”  fry*  + (fi"  h^^)^  y3 

t ~ ( (fi')^  + (x,t)y  + h2(x,t)  (12) 

where  g2  and  h2  are  functions  of  integration.  If  the  solutions  for  (}>i  and  O2 

are  substituted  into  the  boundary  condition,  equation  (6),  where  it  should  be 

noted  that  evaluation  of  any  function  at  y implies  expansion  using  equation 

w 

(7),  then  one  finds  equations  for  g2  and  h,j^.  Thus, 

g2(x,t)  = f h,x  - (-•  ~)  (fi"  (13a) 

± {(2  iz  - (fi')V(V  + 1 ) + (3,(t)  - (f,")^/3}‘^'  (13b) 

where  the  equation  for  h,x  has  been  integrated  once,  and  p(t)  is  the  resulting 
function  of  integration;  f,  and  f2  are  both  functions  of  x and  defined  in  equation 
(7).  It  is  seen  from  equation  (12)  that  in  order  to  complete  the  solutions  to 
second  order,  it  is  necessary  to  find  h2(x,t).  This  is  accomplished  by 


substituting  the  solutions  for  4>i  and  4>2  into  equation  (9c),  integrating  once  to 
obtain  <t>jy,and  substituting  the  resulting  equations  into  the  boundary  condition. 
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equation  (6).  One  finds  that 

r 7 (\  + I ) ]''/360  + (f,'  f/"  - - (f,'  - fi  f," 

+ [(3/2  - - V (f,")^]/3  + f,"  (f," 

- + Pz(t)]/h,x 

where  P^lt)  is , again,  the  function  of  integration.  Finally,  then,  the 
solutions  for  u and  v are,  to  second  order, 

U = 1 + e(f,"y+h,x}4  eM(\+  l)[(f,")"]"  yVZ4  + (\  + l)(f,"h,x)xx 

f ((f,')^+  h,^x)x  + (f,'h,x  - (^*"^>x>x>x  y ^ 

V = $ =e  f,'  + eM(v+ l)[(f,")^]' yV6  + (V+ 1)  (fi"h,x)  yVz 
y X 

+ ( (fl'  + 2 '^'xJx  y " 2 " + • • • 

where  hix  and  h2x  are  given  by  equations  (13b)  and  (14)  respectively,  and 
01  and  02  are  to  be  determined  from  the  initial  conditions  of  the  problem 
under  consideration. 

Before  proceeding  to  the  consideration  of  flows  with  shock  waves, 
it  is  of  interest  to  establish  the  range  of  flow  patterns  covered  by  the 
solutions  represented  by  equations  (15).  First,  it  should  be  noted  that  in 
hix  (equation  (13b)  ),  the  term  0i  - (fi")^/3  cannot  be  negative  since  the 
remaining  terms  under  the  square  root  are  zero  at  x = 0.  If  we  first  con 
sider  the  case  where  0i  - (fi")^/3  > 0 and  also  consider  the  solution  for  u 
(equation  (15a)  ) only  up  to  terms  of  order  €,  then  if  the  upper  ( + ) sign  is 


1 


(14) 


(15a) 


(15b) 


chosen  for  all  x in  equation  (13b),  the  resulting  flow'  is  supersonic  except 
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for  a possible  subsonic  pocket  existing  in  a region  around  x = 0,  as  shown 
in  figure  2a.  If,  on  the  other  hand,  the  lower  (-)  sign  is  chosen  for  all  x 
in  equation  (13b),  one  obtains  a subsonic  flow,  except  for  a possible 
supersonic  pocket  around  x = 0,  shown  in  figure  2b.  In 
both  cases,  the  pocket  enclosed  by  the  sonic  line  grows  as  (3,  - (f,")^/3 
decreases.  Until  Pi  - (fi")^/3  = 0,  no  sign  change  at  x = 0 can  be  tolerated 
because  the  sonic  line  (u  = 1 ) does  not  pass  through  (0,0).  W hen  P i - (f]")^  / ^ = 0 , 
depending  upon  whether  the  sign  in  equation  (13b)  changes  or  remains  un- 
changed as  X passes  through  zero,  the  sonic  line  either  returns  to  the  same 
wall,  figures  3a  and  3b,  or  crosses  the  channel  to  the  opposite  wall,  figures 
4a  and  4l>.  The  izradients  In  velocity,  and  thus  in  pressure,  temperature, 
etc.  , seen  in  these  flow  pictures  illustrate  the  very  interesting  changes 
which  take  place  in  the  flow  when  a radius  of  channel  curvature  of  order  e"’ 
is  considered  as  opposed  to  the  radius  of  curvature  of  order  e considered 
in  reference  6. 

7 

It  can  be  shown  that  for  x = 0(e  ),  the  terms  of  order  e and  in  u 
may  become  of  the  same  order,  casting  doubt  upon  the  validity  of  the 
asymptotic  expansions  in  this  region.  However,  as  shown  in  references 
4 and  5,  this  possible  lack  of  uniform  validity  for  the  asymptotic  expansions 
is  not  a serious  problem,  the  difficulty  disappearing,  for  example,  for  walls 
which  vary  as  x^  in  the  region  where  x = 0(e).  Hence  the  throat  region  will 
not  be  considered  in  this  paper. 

The  solutions  presented  so  far  were  derived  for  flows  with  no  shock 
waves.  It  is  of  interest  now  to  analyze  cases  where  shock  waves  exist,  to 
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ascertain  if  and  at  what  order  these  solutions  break  down.  We  cht^ose  for 
illustrative  purposes  those  flows  for  which  Pi  - (fi")^/3  = 0 shown  in  fi^iures 
3b  and  4a.  In  these  cases,  in  eqxiation  (13b),  the  upper  (i  ) sign  (supersonii 
flow)  is  used  downstream  of  the  throat  but  upstream  of  the  shock,  with  t)ie 
lower  (- ) sign  (subsonic  flows)  holding  downstream  of  the  shock;  the  sign 
upstream  of  the  throat  depends  upon  whether  the  flow  is  subsonic  (accelerating, 
figure  4a)  or  supersonic  (decelerating,  figure  3b).  Similar  results  hold  for 
other  sign  combinations. 

If  the  shock  position  is,  to  lowest  order,  + • ■ • > then  from 

equations  (la)  and  (^),  it  is  seen  that  the  shock  velocity,  dx^/dT,  is  of 
order  . Hence,  to  order  e,  the  shock  polar  equation,  which  reduces  to 
the  Prandtl  relation,  may  be  written  in  terms  of  absolute  velocities.  Thus, 


u j = - u, 

1 d lu 


(lb) 


where  Uj  = 4>ix>  the  sut^scripts  u and  d denote  values  immediately  up- 

stream and  downstream  of  the  shock  wave,  respectively.  Now,  if  equations 
(1  1 ) and  (13  b)  are  used  to  w'rite  U|.^j  and  u^  ^ at  x^ , w'he  re  the  upper  sign  in 
equation  (1  ib)  is  used  for  u^^  and  the  lower  for  u^^j,  it  is  seen  that  the  shock 
w'ave  jump  condition,  equation  (16),  is  not  satisfied.  Thus,  even  in  first 

order,  the  shock  wave  jump  conditions  are  not  siitisfied,  as  opposed  to  the 

6 

case  for  flow's  with  smaller  curvature,  and  an  inner  region  about  the  shock 
must  be  considered. 


SOLUTIONS  IN  THE  SHOCK  WAVE  REGION 
In  the  inner  region  containing  the  shock  wave,  the  solutions  must 
satisfy  the  jump  conditions  across  the  shock,  and  match  with  the  outer 
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channel  flow  solutions  in  the  appropriate  limit.  This  region  may  be  expe<  ted 
to  extend  across  the  channel,  but  to  be  very  thin  in  the  flov  direc  tion.  Then, 
in  the  inner  region,  we  define  the  independent  variables  as  follow's; 

X - X()(t) 


T'“  = T = T f 


f 1 7 a , b , c ) 


where,  again,  Xo(t)  is  the  shock  location  to  lowest  order,  and  cr  is  to  be 
determined.  The  inner  region  velocity  potential  is  defined  as, 

#'  (k'  , y*,  t")  = $(x,y,T)/(V  ^ - .x,,  x^''  (18) 

where  Xg  “ dxg/dT  = ke^dxg/dt.  That  is. 


U - Xf 


1/2,0- 


, v'"  = $'■'  = v/(V  + 1 

y 


( 1 Oa , b I 


In  view  of  the  expansion  for  $,  equation  (4a),  we  write  the  expansion  for  <J>' 
as  follows: 

(x=^  y^  T*)  = (I  -xg)x*  4 (x=^  T*  ) 4 Xg/l'f  4 (20) 

where  ct)''  (x'  , y"' , T’’" ) is  the  inner  perturbation  potential. 

The  proper  form  for  the  expansion  of  4)'  results  from  the  application 
of  three  conditions  which  must  be  fulfilled.  Since  it  is  expected  that  the 
fluid  acceleration  in  the  inner  region  is  of  importance  and  will  be  found  in 
the  lowest  order  governing  equation  for  , in  order  that  changes  intro- 
duced by  the  shock  may  be  accommodated,  the  first  condition  is  that 

A <t>"  * be  of  the  same  order  as  4>?  =:=  =;=  • The  next  condition  is  given 

by  the  matching  conditions  which  are  to  be  met  as  | x"'' | -*oo;  they  are  found 
by  expanding  the  outer  solutions,  equations  (15),  about  x = x^^  and  sub- 
stituting the  resulting  equations  written  in  terms  of  inner  variables,  into 
equations  (19).  I'he  final  condition  is  imposed  by  the  boundary  conditions 


r 


1 1 


at  the  wall.  Since  y - = 0 at  the  walls,  then  also  the  Eulerian  derivative 

of  this  function,  ^{y-yJ/^T  = 0 at  the  walls.  This  equation,  which  may 
be  used  to  derive  equation  (6),  is  transformed  using  equations  (17),  (19), 
(20),  and  equation  (7)  expanded  about  x^,  such  that  it  is  valid  in  the  moving 
coordinate  system  associated  with  the  inner  region,  and  gives  the  desired 
boundary  conditions.  The  three  conditions  lead  to  cr  = 1/2  in  equation  (17a), 
and 


(x*.  y'',  T")  = y*/(Y  + (x*\  y*  . ) 

where  fi;,  = fi'(x(,). 

The  governing  equation  for  in  the  inner  region  may  be  derived 
from  equation  (5)  by  using  equations  (4a),  (18),  (20),  (21),  and  (17)  with 
a - 1 /2 . Thus  , 

*lx*  '*’lx*x*'  ‘*’ty*y*  " 0 


(21) 


(22) 


The  matching  conditions  to  be  used  as  | x*  | -*  oo  , found  as  described  pre- 
viously, are 


f'x* 

ly*  10 


(23a, b) 


where  hj^^  (xq)  is  given  by  equation  (13b)  with  the  upper  (+)  sign  and  where 
the  upper  and  lower  signs  in  equation  (23)  correspond  to  the  conditions  as 
x*-*-oo  andx*-*  + oo  respectively.  The  boundary  conditions  at  the  wall, 
again  found  as  described  previously,  are  as  follows: 

4.,*y*(x^±l,  t=*')  = f,^x*  (24) 

The  final  condition  necessary  before  solutions  in  the  inner  region  may  be 
found  are  the  shock  wave  jump  conditions.  Now,  the  shape  of  the  shock. 
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X t'*'),  is  given  in  the  inner  region  by  the  equation 


dx 


V - V 

d u 

Jic  Jjt 

U J - u 
d u 


(25) 


In  view  of  equation  (21),  it  is  seen  that  the  first  term  in  the  expansion  for 
X*  is  of  order  unity  and  the  next  is  of  order  . In  the  inner  and  outer 
variables,  then,  the  expression 

X (y,  t)  - X 


^ = X *(y*,  t*)  = X,  (y^t")  + (y^t 


* i. 


(Y+ 


.1/2. 


* .* 


) + 


(26) 


where  Xj  (y*,  t*)  = xj  (y,t),  etc.  , meets  the  necessary  conditions  and  agrees 
with  the  relation  used  previously,  that  x^  = Xg  + ....  The  function  x,  is 
chosen  such  that  Xj  (- 1 , t)  = 0 ; that  is , x^  gives  the  position  of  the  inter- 
section of  the  shock  wave  and  the  lower  wall,  and  X(  gives  the  variation  of 
the  shock  shape  from  the  normal  (to  the  x axis)  at  that  point.  In  order  to 
find  the  jump  conditions  across  the  shock  wave,  it  is  necessary  to  write 
the  shock  polar  equation  relative  to  the  moving  shock.  The  velocity  com- 

A A 

ponents  relative  to  the  shock  wave  are  u = u - u , and  v = v - v where 

s n sh 

u , and  V , are  the  components  of  the  shock  wave  velocity.  It  can  be  shown^ 


that  u . = Xq  + 0(6®'^^)  and  v = and  since  x„  = e^kdx  /dt  = O(e^), 


sh 


sh 


it  is  seen  that  if  u,’'  = v,*  = 4>|**,  etc.  then 


ly’ 


u = 1 + e u,*  + u * 


V = (V  + 1 (e  V*  + V,*  + • . . ) 

1/2 


(27a) 

(27b) 


That  is,  to  the  order  desired,  the  shock  wave  velocity  does  not  appear. 


Finally,  from  the  energy  equation  which  holds  in  d moving  coordinate  system. 


it  can  be  shown  that  the  dimensionless  critical  speed  of  sound  is,  in  this 
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system, 

= 1 - 2 (Y  - l)Xo/(V+  1)  + • • • (28) 

The  jump  conditions  are  found,  then,  by  expanding  equations  (27)  about  x*. 
using  equations  (26),  and  substituting  these  expansions  and  equation  (28) 


with  Xq  = e^kdxg/dt,  into  the  shock  polar  equation.  The  results  are  ac 
follows 


It  may  be  noted  that  the  form  of  expansion  derived  for  u*  and  v*  using 


(29) 


boundary  and  matching  conditions  is  also  consistent  with  the  shock  jump 


conditions . 

7 

It  is  not  difficult  to  show  that  upstream  of  the  shock  wave,  the 
solutions  for  are  simply  the  continuation  of  the  outer  solutions  evaluated 
at  the  shock  wave.  That  is,  one  obtains  the  physically  correct  result  that 


no  disturbances  from  the  shock  affect  the  supersonic  flow  upstream  of  the 


shock  wave.  Downstream  of  the  shock,  however,  equation  (22)  must  be 
solved  with  the  boundary  conditions  given  by  the  condition  at  the  walls, 
equation  (24),  the  matching  condition  as  x*-»oo,  equation  (23a),  and  the 
conditions  immediately  downstream  of  the  shock  wave,  given  by  equations 
(25)  and  (29).  Thus,  since  both  u^^  and  v^^  are  found  in  the  first  order 
jump  condition  given  by  equation  (29),  the  shock  wave  is  not  planar  to  this 
order,  and  another  equation,  (25),  is  necessary  to  relate  the  local  shock 
slope  to  Vjj  and  u^^  . It  may  be  noted  that  this  problem  is  therefore  quite 
different  from  those  solved  previously  for  channels  with  smaller  curvature, 
where  the  shock  is  planar  to  the  lowest  order.  In  that  case,  the  boundaries 
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and  normal  derivatives  of  the  potential  at  the  boundaries  were  all  given,  and 
'id  calculated,  as  was  the  higher  order  shock  shape.  In  the  present 

case,  however,  Uj^  and  v^^  are  related  to  the  shock  shape  which  must  be 
found  as  part  of  the  first  order  solution. 

It  is  convenient  to  write  (j)^*  in  terms  of  the  difference  between  it  and 
its  limiting  value  as  x oc  . Thus,  we  define  ?>*as  follows: 


- <t>i*  - tfi'J  y''-  h^^  (Xq)]x 


Then  the  governing  equation  and  boundary  conditions  for  the  problem  to  be 
Solved  are 


^ s.'t  * 

X XX  y y 


where 


* 

X = Xj 


X “*"00 


y*  = ± 1 


1 


r * - ? f " * L ■>  / \ 

- 2f.oy  .2  (— ) 

. 9xi 


, 9x 


i \2 


(3  la) 
(3  1b) 


» * , / PXj  ,r  n A ,9xi.2t 


3y 


3y 


y* 


9x, 


9x, 


(Y^  1) 


9y  %y 


6y 


dy 


^ >,v 

V J - V 

—id _UL 

u * - u * 
id  lu 


(ilii 
(3  ’ 

(32al 


and  equations  (15)  evaluated  at  the  shock  wave,  and  equations  (19),  have  been 
used  to  write  u^*  and  v^*.  Evidently,  this  problem  must  be  solved  numerically 
specific  example  solutions  are  discussed  later. 

A composite  solution,  uniformly  valid  throughout  the  channel,  may  be 
constructed  by  adding  the  outer  and  inner  solutions  and  subtracting  those 


]5 


(matching)  terms  common  to  both.  In  view  of  the  definition  of  t,  , these 
composite  solutions  are  thus,  to  first  order, 

u = 1 t e (fi"  y ± h,„  ♦ ^ + • • • 


V = e fi'  4 e I 1 )‘  ^ 


,1  /z  r * 


(33a) 

(33b) 


where  hj^  is  given  by  equation  (13b)  with  the  plus  sign,  and  = 0 upstream 
of  the  shock,  x < x^.  Also,  for  x < x^,  if  the  flow  accelerates  from  sub- 
sonic to  supersonic  flow  (with  the  sonic  line  passing  through  (0,0)),  then 
the  upper  sign  is  taken  in  equation  (33a);  if  the  flow  is  supersonic  decelerating 
flow  for  X < 0 and  accelerates  for  x ■ 0 with  possibly  a subsonic  pocket  in  the 
region  about  x = 0,  then  the  lower  sign  is  taken  for  x < 0 and  the  upper  for 
X > 0.  For  X ''  Xg,  only  the  lower  sign  in  equation  (33a)  is  taken.  It  should 
be  noted  that  for  x*  = 0(1'  equations  (33)  reduce  to  the  inner  solutions  down- 
stream of  the  shock,  while  as  x * -•  oo  , -*0  and  the  outer  solutions,  to 

order  e , are  recovered. 

Solutions  for  the  pressure  and  density  are  written  in  terms  of  the 
solutions  for  the  velocity  components  through  equations  derived  from  the 
conservation  equations  and  the  equation  for  the  jump  in  entropy  across  the 
shock  wave.  Thus,  from  the  inviscid  energy  equation,  it  can  be  shown  ^ that 


the  stagnation  or  total  enthalpy  can  be  written  as, 

4 h^  4 • • . 


A 

, T I , , p,  V + 1 

\ * 7 ' '' ' = TiFTT) 


= k — (C,%-  2h,^(x„)) 


dt 


ht,  = 0 


X ^ X 


X < x 


(34a ) 
(34b) 


where  h^^  (x^)  is  calculated  from  equation  (13b)  with  the  plus  sign.  The 


(34c) 
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pressure  and  density  are  related  to  the  entropy  and  the  temperature  as 
follows ; 


p/pY  ^ 1 ) 


„ A 

P = p T 


(35a , b) 


The  entropy  change  of  importance  is  that  across  the  shock  wave,  the 


gradients  downstream  of  the  shock  being  of  high  enough  order  that  they  are 
negligible.  For  transonic  flow,  across  a shock 


"d-"  = 


^Y 


2 Y‘ 

- , , , , (M  ^ - 1 f * (M  ^ - 1 )^  + ■ . . 

u 3 (Y  + 1 ) u (Y  + 1 r u 


(36) 


where,  since  the  shock  in  this  case  is  in  motion,  is  the  relative  normal 
Mach  number  of  the  incoming  flow.  For  a shock  described  by  equation  (26), 


0 

1 =€(v*  !)!»,„-  (-5^1*1 


137  ) 


where,  as  mentioned  earlier,  outer  solutions  are  valid  upstream  of  the  shock 
wave.  The  expansions  for  P and  p can  then  be  found  from  equations  (34a) 
and  (35),  where  u and  v are  found  from  equations  (4)  and  (8)  with  u-  = and 

'"i  = *iy 

P = 1 - e Yu,  - €^Y  (v,V2  + U2)  - e^Y{2  (y  4 D [u,^  - 0x,/9y)^]^  /3 

- h^^  - (Y  + 1 ) u,^ / 6 - u,  v,^/2  + V|  v^  + U3}  + • ■ (38a  ) 

p = 1 - eu,  - eM(Y-  l)u,V24  v,V2tU2]  -eM2Y(Y4  1)(U|^  - (dKjbyff/i 

- hj.  j - Y (2  - Y)u,*/3  + (Y  - 1 ) u,  U2  - (2  - Y)  u,  v,^/2 

+ V,  V2  + U3]  + • • • (38b ) 

If  the  composite  solutions  are  used  for  u,  and  v,  (terms  of  order  e in  equations 
(33)  ),  then  equations  (38)  give  solutions  which  are  uniformly  valid  to  order  e . 

A solution  for  ^ similar  to  those  given  for  P and  p is  easily  derived  from 


equation  (34a). 
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It  should  be  noted  that,  if  composite  solutions  are  considered,  then 
immediately  upstream  of  the  shock  wave,  = - e "Y  fjJ,  while  downstream 

Thus,  downstream  of  the  wave. 


of  the  wave,  P, 


if  the  wave  is  approached  by  first  going  to  y*  = ± 1 > and  then  going  to  x * = x,  , 
then  at  the  wave  P^  = -eV  f^'J  . However,  if  the  wave  is  approached  by  going 
to  x*=  xj  first,  at  any  y*,  and  then  traveling  along  the  wave  toward  y*  = ± 1 , 
Py  = t ^ + ■ ■ ■ • 1 hat  is , from  equations  (31c)  and  (3  1 e ) , 9xi  / 3y  = 0 at 

* » t 

y = i.  1 ; hence,  from  equation  (3  lb),  * = - 2f.o  at  y = i.  1 . Thus,  there 

x'V 

is  a change  of  sign  in  P across  the  shock,  and  since  P„  oc  v ’*  *,  the  stream- 

r-  r y y 1 X 

line  curvature,  this  means  that  downstream  of  the  shock  wave,  the  limiting 
value  of  V * * at  the  intersection  of  the  shock  wave  and  the  walls  depends  on 
the  path  of  approach  to  the  intersection.  This  singularity  in  the  solution  is 
similar  to  that  analyzed  by  Messiter  and  Adamson. 


SHOCK  WAVE  LOCATION 

The  unsteady  perturbations  in  flow  velocity  and  thermodynamic 
properties  in  the  channel  can  be  caused,  for  example,  by  pressure  distur- 
bances in  a plenum  downstream  of  the  shock  wave.  Here,  we  consider  the 
case  where  the  disturbances  have  a characteristic  time  (inverse  of 
characteristic  frequency)  of  order  and  an  amplitude  of  order  c^.  As 

will  be  seen,  this  is  sufficient  to  give  shock  motions  with  an  amplitude  of 
order  unity. 

Since  only  terms  of  order  in  P may  vary  with  time,  it  is  seen 
from  equations  (38a)  that  Ui  is  independent  of  time;  hence,  h|  = h,  (x).  In 
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fact,  (equations  (38a)  and  (15))  the  time  dependence  first  appears  in  U2  , 
through  = h2(x.t).  As  a result,  in  equation  (13b),  p,  is  a constant,  and 
the  time  dependence  of  h2  and  thus  U2  and  indeed  ^2  is  through  P2  (i)-  Both 
Pi  and  P2  may  have  different  values  upstream  and  downstream  of  the  shock 
wave.  Since  the  time  enters  the  solutions  for  u only  through  the  integratif)n 
function,  P2,  then  to  order  , the  unsteady  motion  may  be  pictured  as  a 
sequence  of  steady  state  solutions  for  u,  each  with  different  downstream 
conditions.  The  velocity  component  v,  to  order  , is  totally  time  indenen- 
dent,  (equation  (15b)).  From  equations  (38a),  (15),  and  (M),  it  is  seen 
that  specifying  the  second  order  variation  in  pressure  at  a downstream 
plenum  location,  is  equivalent  to  setting  P2  (t)  downstream  of  the  shock, 
say  The  fact  that  time  does  not  appear  explicitly  in  the  solutions  valid 

downstream  of  the  shock  wave  (e.g.  , see  the  composite  solutions  given  in 
equations  (33)  ) and  appears  only  as  a parameter,  in  P2  , is  due  to  the  fact 
that  the  disturbances  travel  upstream  at  a speed  of  order  e while  the  shock 
moves  at  a velocity  of  order  6^  , and  indeed,  the  order  of  the  unsteady  part 
of  the  fluid  velocity  is  . Thus,  the  disturbances  travel  upstream 
"instantaneously"  compared  with  the  characteristic  times  under  study. 

It  should  be  noted,  however,  that  as  a result  of  the  motion  of  the  shock 
(i.e.,  since  x^  = first  order  pressure,  velocity,  etc.  , in  the  range 

of  motion  of  the  shock,  do  change  as  the  shock  moves  back  and  forth;  they 
jump  between  the  steady  state  values  upstream  and  downstream  of  a shock 
wave  at  the  point  in  question. 

The  equation  from  which  the  instantaneous  shock  wave  location  can 
be  calculated  is  derived  from  the  principle  of  mass  conservation  applied 

1 


L 


to  a control  volume  enclosing  the  shock  wave.  Thus,  one  surface  of  the 


control  volume,  the  sides  of  which  are  the  channel  walls,  is  located  at  x - 0, 

the  throat,  and  the  other  is  at  x = X , somewhere  downstream  of  the  shock. 

c 


Then 


X 


d(r^Wc*’c  ^ r r’  vj  ''  c' 

— j pdxdy>-j  p udy  i J p udy  = 0 

/ J /V  \ 


w c 


(39) 


y'  (X  ) 
w c 


where  y and  y denote,  respectively,  the  upper  and  lower  walls.  The 
w w 

integral  in  the  first  bracket,  from  zero  to  X , is  evaluated  in  two  parts, 
one  upstream  of  the  wave  (0  to  x^)  and  one  downstream  of  the  wave  (x^  to  X^) 
with  the  composite  solutions  used  for  U[  in  the  first  order  terms  of  equation 
(38b).  Since  d/dT  = O(e^)  and  only  terms  up  to  third  order  are  desired, 
only  the  first  order  term  in  p is  necessary.  In  the  remaining  integrals  at 
x = 0andx  = X , equations  (15)  and  (38b)  are  used.  First  order  term.s  are 
identically  zero,  and  from  the  second  and  third  order  terms,  one  finds  that 


P„.  =f>, 


(\  + I ) dt 


2d"  10 


10  ' 


10 


9x 


9X|  ^ 1 8xi 


Y r 2 r •>  1 '^'^1  1 


(40a) 

(40hl 

(40c) 


■10  1 "iu'  -o''3y'  3 ' 3y 

where  hjg  = dh,/dx  evaluated  at  x^and  is  a constant, while  P^^  is  a function 
of  time  for  the  case  under  consideration  where  disturbances  are  imposed 
downstream  of  the  shock  wave.  If  the  shock  wave  has  a higher  order  curva- 
ture, 9X|/3y  = 0,  and  G (x„)  vanishes;  then  equation  (40b)  is  reduced  to  a form 


equivalent  to  that  given  by  Richey  and  Adamson. 
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Equation  (40b)  may  be  used  in  determininjj  the  shock  location, 

Xg  = constant,  for  a flow  which  is  steady  with  a stationary  shock  wave.  Thus, 
for  a given  (now  constant)  which  is  equivalent  to  a given  pressure  in  a 
plenum  downstream  of  the  shock  wave,  the  only  unknowns  in  equation  (40b  I, 
h 10  and  G (xg ) , are  functions  ofxg.  Since  8x[/9y,  which  depends  on  Xg  as  veil 
as  y as  seen  from  equations  (31),  must  be  known  if  G (x^  ) is  to  be  calculated, 
this  means  that  equations  (31)  and  (40)  must  be  solved  simultaneously.  In 
principle,  therefore,  x^  may  be  calculated  for  a given  P^^j  . However,  it  is 
much  easier  to  choose  a given  Xg , solve  equations  (31)  and  in  the  process 
find  X(  , calculate  G (x^)  and  hj^,  and  then  use  equation  (40b)  to  calculate  the 
^2d  corresponds  to  the  chosen  Xg ; after  a series  of  such  calculations, 

x^  may  be  plotted  as  a function  of  p2(j-  the  downstream  pressure  varies 
with  time  in  a pr  sr  -ibed  manner,  then  p^^j  = known  and  dx^/dt  and 

thus  the  shock  velocity  may  in  principle  be  calculated  from  equations  (31) 
and  (40b),  v'ith  the  shock  position,  x^(t),  being  found  then  by  integration  of 
dx^/dt.  The  computations  may  be  carried  out  by  first  solving  equations  (31) 
for  a sequence  of  values  for  x^,  so  that  0X|/9y,  h^'^  and  thus  G (x^)  are 
essentially  known  as  functions  of  Xg.  With  p^^(t)  and  the  initial  value  of  Xg 
known,  equation  (40b)  may  then  be  integrated  numerically,  thus  relating  Xg 
to  t;  i.e.  , Xp(t)  is  obtained . Finally , with  x^j(t ) known , ?,*and9xi/9y,  and 
thus  u,  P,  p,  Xi  , etc.  , may  be  obtained  as  functions  of  time  at  any  space  point. 

NUMERICAL  CALCULATIONS 

In  the  present  work,  equations  (31)  were  solved  using  the  method  of 

9 1 0 

integral  relations  proposed  by  Dorodnitsyn.  ’ Details  of  the  computation 

i 

L.  - — d 


Z1 


are  avai 


lable  in  reference  7,  so  only  a brief  description  is  given  here.  The 


inner  region  in  which  equations  (3  1)  are  valid  is  composed  of  a region  extend- 
ing from  y*  = - 1 to  y*  = + 1 in  the  y*  direction  and  from  the  curve  x ’*'=  to 

X*  - 00  in  the  x*  direction;  x = x,  + • • • is  to  be  found  as  part  of  the  solution. 

s ^ 

. A _ -X* 

This  region  is  transformed  into  a finite  region  by  the  transformation  x - e , 

A A X 

and  this  finite  region,  between  x = 0 and  x = x^  = e , is  divided  into  N 

strips.  Calculations  of  the  shock  shape  made  with  N = 2 and  compared  with 

the  calculation  made  with  N = 3 showed  very  little  difference,  indicating  that 

relatively  accurate  computations  can  be  made  with  two  strips.  Equation  (31a), 

written  in  terms  of  velocity  components,  and  the  corresponding  ir rotationality 

condition,  u**  = v**,  are  the  governing  equations  integrated  across  each 
ly'*'  ix^ 

strip,  using  order  polynomial-s  in  x/x^  for  interpolation  expressions  for 

lim  ij, 

Ui*  and  V*.  In  the  present  case,  it  was  found  that  if  = 0 is  used 

as  a boundary  condition  and  the  method  of  integral  relations  applied,  then 

4 **  is  not  zero  as  X*  — or.  That  is,  apparently  due  to  the  approximations 

inherent  in  the  method  and  the  effects  of  the  singularity  at  the  shock-wall 
interaction,  there  is  an  error  in  the  v velocity  as  x*  — oo  . Rather  than  accept 
this  error,  the  condition  as  x*  -*oc  w'as  enforced  and  the  ir  rotationality 

condition  was  integrated  once  across  the  whole  region  rather  than  across 
each  of  the  two  strips  separately. 

Using  the  two  strip  method,  one  finds  a nonlinear  differential 
equation  for  the  shock  shape;  since  time  enters  only  as  a parameter, 
this  equation  may  be  treated  as  an  ordinary  differentia]  equation.  Thus,  if 


7,(y)  = 9xi/9y,  one  finds  that 


Z2 


(fio'y  ^ Vo  - 3z^)z"  - [6zz'  - 2 Rz  + f/o'hj'o  (1  - /)  - 2 f,”]  z> 

- [2R^  ^ (3h,'o  * 2f/Jy)R  - Zf/Jh/gyjz 

- ^ fi'o'  h,'o  (1  - y^)  + 3 f ;■  ] R - 2 f,'J  y - f/J  . 0 

R = [(fi’i)'yV  f;;h;„y  4 (h,-/ . -y^)z  - f.-z 

- (f,”y  ^ hi,  - 3z^)z']‘/^ 

where  h,Q  = hj'(Xy),  etc.  The  boundary  conditions  for  equations  (4‘4)  are 

z (±  1 ) - 0 

Equation  (41a)  uas  solved  numerically  using  a fourth-order  Runge-Kutta 
method,  yielding  z(y).  Then  z(y)  - 9x,/9y  was  integrated,  using  the  con- 
dition that  Xi(-1  ) - 0,  to  give  the  shock  shape,  Xj  . Since  both  f['g  and  h,'g 
are  time  dependent,  because  Xg  = Xg(t),  the  coefficients  in  equations  (41) 
have  different  values  at  different  times  and  it  is  seen  that  the  shock  shape 
varies  with  time. 

The  wall  shapes  chosen  for  the  calculations  are  as  follows; 
y = ± \ 4 e xV2  ± eMl  * 20  (Y  f 1 )/9]  xV2 

w 

where  the  wall  shape  functions,  f[(x)  and  f^(x)  are  found  by  comparing 
equations  (7)  and  (43).  In  order  to  make  the  calculations  as  simple  as 
possible,  Xg(t)  is  prescribed  here  and  the  coi  responding  (and  thus 
pressure  in  the  downstream  plenum)  is  calculated  from  equations  (40b) 
and  (40c).  That  is,  one  can  prescrilse  P^d^*"^  find  the  resulting  Xg(t), 
or  prescribe  Xg(t)  and  find  the  necessary  p^^(t).  The  latter  problem 
demands  less  coniputing  time  and  is  used  here  for  illustration.  Thus, 

X (t)  = 1.54  (sin  t)/2 
0 


(41a) 


(41b) 


(42) 


(43) 


we  set 


(44) 
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and  perform  calculations  for  V = 1 . , e “ 0.  1,  and  = 0.  The  wall  shapes 

given  by  equations  (43)  with  these  values  for  e and  V,  are  those  seen  in 
figures  2 to  4. 

Velocity  profiles  showing  the  first  order  inner  solutions  downstream 
of  the  shock  wave  are  shown  at  x * = x,  in  figure  5a , and  at  x * = X(  + 0.  69  in 

TT 

figure  5b,  each  at  three  different  times , ^ ~ ^ . The  corresponding 

is  found  from  equation  (44).  These  figures  thus  illustrate  the  temporal 
variations  in  the  velocity  components  at  the  indicated  inner  region  stations. 
The  relatively  large  spatial  variations  which  occur  in  the  inner  region  at  a 
given  time,  are  illustrated  in  figure  6 , at  t = 0 . The  change  in  Ui*  from  its 
values  immediately  behind  the  shock  to  the  linear  profile  associated  with 
subsonic  flow  for  the  given  wall  shapes,  is  indicative  of  the  large  accelera- 
tions and  decelerations  which  take  place  in  the  inner  region. 

The  shock  shape , x^  = Xg  + e ' (V  + 1 Xj  + • • • , is  show-n  in  figure  T , 
as  a function  of  time.  It  is  seen  that  as  a result  of  the  curvature  of  the  walls 
and  the  attendant  gradients  in  the  incoming  flow,  the  shock  w'ave  has  a pro- 
nounced curvature;  for  the  direction  of  curvature  chosen  for  the  walls,  the 
shock  starts  normal  to  the  lower  wall,  inclines  in  the  flow  direction  and  then 
turns  back  toward  the  upper  wall  so  as  to  become  normal  to  it  also.  As  the 
shock  moves  closer  to  the  throat  of  the  channel,  it  becomes  weaker  and  has 
more  curvature. 

The  values  of  Xg,  for  each  time  at  which  calculations  were  made,  the 
cor repsonding  calculated  from  equation  (40b),  and  the  equivalent  change 

in  back  pressure  at  a downstream  plenum  from  its  value  at  t = 0 , ~ 


Pb^®^>  Siven  in  table  1.  Finally,  the  pressure  on  the  upper  and  lower 
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walls  at  these  same  times,  to  order  e . are  shown  in  figure  8.  In  this  figure, 
the  solid  lines  indicate  the  calculated  pressure  distributions  upstream  of  the 
shock  wave,  the  pressure  jump  across  the  shock,  and  outer  pressure  distri- 
butions downstream  of  the  shock.  The  large  dots  indicate  the  pressures 
calculated  in  the  inner  region,  and  the  dashed  lines  show  a curve  drawn 
through  these  points  and  faired  into  the  outer  pressure  distribution,  to  shou 
an  approximate  pressure  distribution.  Since  only  two  strips  vere  used  in  the 
inner  region  computations,  only  one  data  point  within  the  inner  region  is 
available . 


CONCLUSIONS 

The  methods  used  in  this  study  enable  one  to  study  unsteady  flows 
with  shock  waves  in  relatively  highly  curved  asymmetric  channels  with 
arbitrary  wall  shapes  and  impressed  disturbances  of  arbitrary  form.  It  is 
shown  that  the  shock  wave  is  not  planar  and  that  its  shape  must  be  obtained 
as  part  of  an  inner  solution  which  involves  a numerical  solution  of  the  non- 
linear small  disturbance  transonic  equation,  the  unknown  shock  shape  form- 
ing one  of  the  boundaries  of  the  region  in  question.  It  proves  to  be  relatively 
easy  to  use  the  method  of  integral  relations  to  obtain  approximate  but  very 
useful  solutions.  The  solutions  allow  one  to  calculate  the  shock  shape  and 
velocity,  as  well  as  fluid  velocity  and  thermodynamic  property  distributions. 


as  functions  of  time. 
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Figure  1.  Sketch  of  asymmetric  channel  flow,  showing 
coordinate  system. 


-i.O 


/ 1 1 / 1 n / / / f / ' • ' ' ■ 


(a)  + sign  for  x < 0 and  - sign  for  x > 0 in  equation  ( I 3b ) 


(b)  - sign  for  x < 0 and  sign  for  x > 0 in  equation  (I3b) 


Figure  3.  Isotachs  corresponding  to  equation  (I  5a)  when 
Pi  = (fi")^/3  = 1/3,  for  same  wall  shapes  and  e 
as  in  figure  2 . 


Figure  5.  First  order  velocity  component  profiles  downs 
times,  for  channel  walls  given  by  equation  (43 
(44).  V = 1 • ‘i- 


component  pr 
hannel  walls  p 


Figure  A.  First  order  velocity  component  profiles  at  various  positions  downstream 
of  the  shock  at  t = 0,  for  channel  walls  given  by  equation  (43)  and  shock 
location  as  in  equation  (44).  o *(or , y*,  0)  = 0.  y = 1.4. 


Figure  H.  First  order  wall  pressure  distributions  for  various  times  for  i hannel 
walls  given  by  equation  (43)  anrl  shock  location  as  in  equation  (44). 

Y = 1 . 4,  e = 0.  1 . calculated  pressure  distribution  upstream  of  the 

shock  wave,  jump  across  shock,  and  outer  pressure  distribution  down- 
stream of  shock.  • pressures  calculated  in  inner  region.  apnroxi- 

mate  pressure  distribution  in  inner  region. 
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